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Abstract — This paper is concerned with the study 
of the peristaltic flow of Herschel - Bulkley fluid in an 
inclined flexible channel lined with porous material under 
long wave length and low Reynolds number assumptions. 
This model may be applicable to describe blood flow in 
the sense that erythrocytes region and the plasma regions 
may be described as plug flow and non-plug flow regions. 
The effect of yield stress , Darcy number , angle of 
inclination and the index on the flow characteristics is 
discussed through graphs. 
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X. Introduction 
Peristaltic transport is a form of fluid transport 
which occurs in biological systems. This mechanism has 
received considerable attention in recent times in 
engineering as well as in medicine. It plays an 
indispensable role in transporting many physiological 
fluids in the body. Many modern mechanical devices 
have been designed on the principle of peristaltic 
pumping for transporting noxious fluids without 
contaminating the internal parts. Further the blood 
transfusion process in dialysis, the mechanism of 
peristalsis may be applicable since the blood flow in 
small blood vessels is reported to be done by peristalsis. 



Latham (1966) made first experimental study of 
the mechanics of peristaltic transport. The results of the 
experiments were found to be in good agreement with 



the theoretical results of Shapiro (1967). Based 
on this experimental work. Burns and Parkes (1967) 
studied the peristaltic motion of a viscous fluid through a 
pipe and a channel by considering sinusoidal variations 
at the walls. 

In physiological peristalsis, the pumping fluid 
may be a Newtonian or non - Newtonian fluid. Kapur 
(1985) suggested several mathematical models for 
pumping physiological fluids. Among these some 
models deal with Newtonian fluids and others with non - 
Newtonian fluids. Ravi Kumar et al. [ ] studied the 
peristaltic pumping in a finite length tube with permeable 
wall. Krishna Kumari et al. [ ] studied the peristaltic 
pumping of by considering Jeffrey model. 

Scott Blair (1959) reported that blood obeys the 
Casson model only for moderate shear rate flows. He 
also reported that the assumptions included in Casson's 
equation are not suitable for Cow's blood and that the 
Herschel - Bulkley equation represents fairly closely 
what is occurring in the blood. Herschel - Bulkley fluid is 
a semi solid rather than an actual fluid. 

Among models of semisolids, the Herschel - 
Bulkley model is preferred because it describes blood 
behaviour very closely and also the Newtonian, 
Bingham and Power - law models can be derived as 
special cases. Furthermore, Herschel - Bulkley fluids 
describe flows with a non -linear stress strain 
relationship either as a shear thickening fluids or shear 
thinning one. Since shear thinning and shear thickening 
fluids play an important role in biomedical engineering. 
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Some examples of fluids behaving in this 
manner include food products, pharmaceutical products, 
slurries, polymeric solutions and semisolid materials. 
Chaturani and Samy (1985) discussed the blood flow 
through a stenosed artery by considering blood as a 
Herschel - Bulkley fluid. The gastrointestinal tract is 
surrounded by a number of heavily innervated smooth 
muscle layers, contraction of these muscle layers can 
mix the contents of the tract and move food in a 
controlled manner in an appropriate direction. Epithelial 
cells, beneath these layers are responsible for the 
absorption of nutrients and water from the intestine. 
These layers consist of many folds and there are pores 
through the tight junctions of them. So the flow of fluids 
in different geometries in channels/tubes with porous 
material at the boundary is very significant in 
physiological applications. Vajravelu et al.[ ] considered 
the Herschel - Bulkley fluid in their study of peristaltic 
pumping of fluids. 

In view of these, peristaltic flow of Herschel - 
Bulkley fluid in an inclined flexible channel lined with 
porous material is studied under long wave length and 
low Reynolds number assumptions. This model may be 
applicable to describe blood flow in the sense that 
erythrocytes region and the plasma regions may be 
described as plug flow and non-plug flow regions. 

HERSCHEL - BULKLEY MODEL 

The basic equations governing the flow of an 
incompressible Navier-Stokes fluid are the field 
equations 

divV = 0, 



diva pf - 



dv 
dr 



where V is the velocity, f the body force per unit mass, n 
the density, and d/dt the material time derivative. 6 is 
the Cauchy stress defined by 

a = -pI-T 

T = 2juD + S, S = 27]D 



2^ 



, L = grad V 



Also -pi denotes the indeterminate parof the stress 
due to the constraint of incompressibility, jU and c are 
viscosities. 

The Herscel-Bulkely model combines the effects of 
Bingham and power-law behaviour in a fluid. For low 

strain rates (j^< Tq ///q), the "rigid" material acts like a 
very viscous fluid with viscosity //^ . As the strain rate 

increases and the yield stress threshold, Tq is passed, 
the fluid behaviour is described by a power law 



7 



where k is the consistency factor and n is the power-law 
index. 



II. mathematical formulation 

Consider the peristaltic pumping of a Herschel - 
Bulkley fluid in an inclined flexible channel of width 'a'. 
The channel is inclined at an angle 9 with the horizontal. 
The flexible walls of the channel are coated with non 
erodible porous material. Suppose we have a 
longitudinal train of progressive sinusoidal waves on the 
upper and lower walls of the channel. For simplicity, we 
restrict our discussion to the half width of the channel as 
shown in Fig. 1. The region between y=0 and y = yo is 
called a plug flow region. In the plug flow region, 

I Txy| ^ . In the region between y = yo and y = H, we 
have I Txy I > To . 



The wall deformation is given by 
Y = H(X, t) = a + b sin (X - ct) 



(1) 



where D is the symmetric part of the velocity gradient, where b is the amplitude , X is the wavelength and c is 
that is, the wave speed. 
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Under the assumptions that the channel length 
is an integral multiple of the wavelength X and the 
pressure difference across the ends of the channel is a 
constant, the flow becomes steady in the wave frame 
(x,y) moving with the velocity c away from the fixed ( 
laboratory ) frame (X,Y) . 

The transformation between these two frames is given 
by 

X = X-ct , y = Y , u(x,y) = U (X-ct, Y) - c 



v(x,y) = V(X-ct , Y) , p(x) = p(X,t) , 0 = ir-Y 



(2) 



— l^vv ) = — ^ + 77 sin 6' 

'dy^ '''' dx 



(3) 



where T = 



V 



dy 



0 



(4) 



J 



The dimensionless boundary conditions are 



y/ = 0 



a{y = 0 



(2) 



(5) 



where U and V are velocity components in the 
laboratory frame and u, v are velocity components in the 
wave frame and © and ir are the stream functions in 
the wave and laboratory frames respectively. In many 
physiological situations it is proved experimentally that 
the Reynolds number of the flow is very small. We 
assume that the flow is inertia free and the wavelength 
is infinite. 

We introduce the following non-dimensional 
quantities in order to make the basic equations and 
boundary conditions dimensionless. 



_ X _ y — h 
X =—; y = —; h =—; 
X 



a 



a ( 



a 



; ¥ = 



- Ct _ 8 

A a 



ac 



ac 



_ u dy/^^^ _ vX dy/ _ pa 
u- — -— — v= — = —\ p- ; 

c d y/y Q ac dx jiXc 

where Ji and V are velocity components in the wave 
frame. 



y/ =0 



at j; = 0 



(6) 



T =0 



aty = 0 



(7) 



a 'dy 



(8) 



III. SOLUTION OF THE PROBLEM 
On solving equations (3) and (4) together with boundary 

dy/ dy/ 

conditions (5) - (8) and u = and V = , we get 



dy 



dx 



the velocity field as 



u = F 



K + 



1 f/, x/c+i / ^ Da / x/ 

-\[h-^-y,) -[y-y,) ] +—[y-y,) 



a 



-1 



where p = vrj^md: k = — 

dx n 



(9) 



(10) 



Under the lubrication approach, the equations governing 
the motion become (dropping bars) 



We find the upper limit of the plug flow region using the 
boundary conditions that y/^^ = 0 at = so, we 



^0 

have, J^o = — 



(11) 
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Also by using by the condition 

^yx ~ ^h-& y = h-e ,\Ne obtain 



h- 



(12) 



Hence = ^ = r;0<r<l 



h-e T 



h- 



Taking y = yQ in (9) we get the velocity in the plug flow 
region as 



" k+l 



(13) 



Integrating equations (9) and (13) and using the 
conditions y^„=0 and ^ = 5^ at j = we get the 

stream function as 



1 



¥=-y+p''{-7^Ah-€-yof"'y 

k+l 



iy-yo) 



k+2 



k+2 



(14) 



¥p=-y+-^\ih-^-yo) 



y 



(15) 



The volume flux 'q' through each cross section in the 
wave frame is given by 



yo 



h-t 



q = \Updy+ \udy 
0 



P 



= -{h-e) + i—[ih-e-y,r\h-e)- 



k + \ 



h-e-y^ ^Da 



k + 2 



+ 



(16) 



a 



The Instantaneous volume flow rate in the 

laboratory frame between the central line and the wall is 



H 



Q{X,t)= \u{X,Y,t)dY 



(17) 



0 



From equation (16) we can write 



a k+\ 



dp 
dx 



{q + h-e){k + \){k + 2)a 



k 



{h-sf' {l-rf' {k + 2)a -{h-Gf {l-rf \{h-G){l-T) + s[D^{k + 2) 



+Tjsin0 



(18) 

Averaging equation (17) over one period yields the time mean flow (time - averaged volume flow rate ) Q as 



^ 0 



(19) 



The pumping characteristics 
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Integrating the equation (18) witin respect to V over one waveiengtin, we get tine pressure rise (drop) over one 
cycle of the wave as 



i 

-J 



{Q-l-{-h-e){k-\-l){k-h2)a 



[{h-ef'\l-Tf\k^2)ay{h-ef\l-Tf'[{h-e){l-T)^^fD^{k^2) 



dx 



(20) 

The tinne averaged flux at zero pressure rise is denoted by 2^ and the pressure rise required to produce zero 
average flow rate is denoted by AP. . 



1 

-I 



{-\ + h-e){k + \){k + 2)a 



[{h-ef\\-T)''\k + 2)ay[h-e)'"\\-4'' [[h-e)[\-T) + 4Da{k + 2) 



dx 



(21) 



The dinnensionless frictional force F at the wall across one wavelength is given by 



dx 



dx 



J 



F = 



1 



h{-l + h-e){k-^l){k + 2)a 



[{h- e (1 - {k + 2)a]-{h-e f"' (1 - [{h-e){l-T) + ^{k + 2) 



+ 77 sin 6^ 



dx 



(22) 



IV. DISCUSSION OF RESULTS 

The variation of pressure rise with time averaged flux is 
calculated from equation (20) and is shown in Fig. (2) for 
different values of yield stress with 

7] = 3,Da = O.Ola = O.OlO = 0,(/> = 0.6,g = 0. 1 and 

77 = 0.2. The pumping curves intersect at a point (0.5, 
0.006) in the first quadrant due to variation in the 
parameters. It is found that AP increases with the 
increase in yield stress Tto the left of the point of 



intersection and opposite behaviour is observed to the 
right of this point. For a fixed AP , the flux increases due 
to an increase in the yield stress parameter r in the first 
half of the channel. 

The variation of APwith Q for different values 

of r for 0 = TT /3 is shown in Fig. (3). We find that the 
behaviour is same for horizontal and inclined channels. 
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The variation of pressure rise with time 
averaged flex is numerically evaluated for different 
values of Darcy number 'Da' with 

n = 3,a = 0.0l,(/> = 0.6,z = 0.1 and is shown in Fig. 
(4). It is observed that for a given Ap , the flux decreases 

with an increase in the Darcy number for 0<Q< 0.5 
and the opposite behaviour is observed forO < 2 < 0-5 ■ 
For a given flux Q > 0.5 , the pressure difference 
AP decreases with an increase in Da and the behaviour 
is found to be opposite for Q > 0.5 . 

The variation of pressure rise Ap with the flux 

Q is calculated for different values of 'g ' with 0 = 7r/3 , 

(/> = 0.6, T = 0.1 and a = 0.01 and is shown in Fig. 
(5.) It is found that for a given 'g ' (thickness of the 
porous lining), Ap decreases with the increase in the 

time averaged flux Q. For a given Q >0.5 
(approximately), Ap increases with the increasing g 

and Q > 0.5 the variation of Ap is negligible. 

The variation of pressure rise with 2 for 

different values of angle of inclination '0' for different 
values of angel of inclination is shown in Fig (6). It is 

observed that for a given Ap , the flux Q increases 
with the angle of inclination of the channel with the 
horizontal. For a given 2 , it is found that the pressure 
rise increases with the thickness of the porous lining. 

The variation of pressure rise with time averaged 
flux is calculated for different values of n (index) and is 



shown in Fig. (7). We observe that the pumping curves 
intersect at a point (0.5, 0.17) in the first Quadrant. For a 

given Ap, we find that the flux Q decreases with an 

increase in index, 'n' for 0 < 2 < 0.5 heng > 0.5 , the 

opposite behaviour is noticed. For a given flux Q < 0.5 . 
The pressure rise decreases with increasing 'n'. The 
opposite behaviour is observed for Q > 0.5 . 



1^ 




~ >x. 




Fig. 1 . Physical model 
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K 



Fig 2. Variation of tSp with Q for different values of T when n 
3,Da = 0.01, 6^ = 0.01, 6> = 0,G = 0.1, ^ = 0.6, 77 = 0.2 
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Fig 3. Variation of Ap with Q for different values of T when 
n = 3,Da = 0.01 , a = 0.01, 0 = 7r/3,e = 0.1,^^ = 0.6,7] = 0.2 ■ 
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Fig 6. Variation of Ap with Q for different values of 0 when 
n = 3, Da = 0.1,6kr = 0.1, G = 0.1, = 0.6,// = 0.1,T = 0.1 
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Fig 7. Variation of ISp with Q for different values of n winen 
Da= 0.2,cr = 0.1,G = 0.1, 6> = ;r/3, ^ = 0.6, 77 = 0.1,7 = 0.2 
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